Understanding spin glass transition as a dynamic phenomenon 
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Existing theories explain spin glass transition in terms of a phase transition and order parameters, 
and assume the existence of a distinct spin glass phase. In addition to problems related to clarifying 
the nature of this phase, the common challenge is to explain profound dynamic effects. Here, we 
propose that the main experimental results of spin glass transition can be understood in an entirely 
dynamic picture, without a reference to a distinct spin glass phase, phase transition and order 
parameters. In this theory, the susceptibility cusp at the glass transition temperature is due to 
the dynamic crossover between the high-temperature relaxational and low-temperature spin wave 
regime. The crossover takes place when t — r, where t is observation time and r is relaxation 
time. Time-dependent effects, inconsistent with the phase transition approach, and the logarithmic 
increase of T g with field frequency in particular, originate as the immediate consequence of the 
proposed picture. In our discussion, we explore similarities between the spin and structural glass 
transitions. 



INTRODUCTION 

Since its discovery, spin glass has been considered as a 
third distinct type of low-temperature magnetic arrange- 
ment in solids, in addition to ferromagnetic and anti- 
ferromagnetic [lj. In spin glass, spins are disordered, 
prompting the term "glass" , and give zero net magnetic 
moment. Great theoretical effort went into understand- 
ing the nature of the spin glass transition. The collection 
of new ideas and theories has formed a new large research 
field in condensed matter physics, with connections pro- 
posed to other disciplines such as economics and biology 
lHH]. A couple of decades ago, it was remarked that "as 
fame and topicality diminish, the problem remains an in- 
teresting, yet passe, research topic" One of the rea- 
sons for continued interest in this area is the presence of 
fundamental theoretical problems that are not resolved. 
There is also current interest in new systems such as mag- 
netic nanoparticles, ferrofluids, relaxor ferroelectrics and 
so on, where spin glass behavior is observed. 

The main signature of the spin glass transition is the 
susceptibility cusp, which is considered to mark the tran- 
sition from the paramagnetic to the spin glass phase [lj]. 
The temperature at which the cusp develops is called spin 
glass transition temperature, T g (or Tf). The suscepti- 
bility cusp has stimulated approaches to spin glass tran- 
sition based on the existence of a phase transition at T g , 
between the high-temperature paramagnetic phase and 
a novel low-temperature "spin glass" phase. This was 
primarily inspired by the phase transition theory, cen- 
tered at the changes of system properties at the phase 
transition temperature 0, @ • 

The traditional phase transition theory requires the 
presence of two distinct phases on both sides of the phase 
transition 0, [j| . On the other hand, as detailed stud- 
ies revealed, there is no obvious distinct second phase 
in the spin glass transition: the structure of the low- 
temperature spin-glass phase is not different from the 



high-temperature paramagnetic phase, and is equally dis- 
ordered. The absence of a second phase presents a prob- 
lem with viewing spin glass transition as a phase tran- 
sition, resulting in persisting difficulties and controver- 
sies [l|. To circumvent the problem, popular theories 
introduced and discussed novel non-conventional order 
parameters [H-Q while retaining the view that a phase 
transition of some sort takes place. An advantage of 
this approach is that well-established methods of statis- 
tical mechanics can be used to calculate system equilib- 
rium properties and their derivatives such as suscepti- 
bility. Such was the Edwards- Anderson (EA) theory 
which introduced the new order parameter for spin glass: 
q = (s\ 8*- ), where superscripts 1 and 2 denote two 
different moments of time for spin i. 

As introduced, q serves to quantify the difference be- 
tween the mobile spins in the paramagnetic phase and 
the frozen spins in the glass phase at T g . Notably, q 
describes the change in dynamic behaviour but not the 
difference between two distinct equilibrium phases as in 
the conventional phase transition theory. The important 
next step in the EA theory is the substitution of the time 
average by the ensemble average, by making an assump- 
tion that the system is in equilibrium and is ergodic. 
Subsequently, statistical mechanics is used to calculate 
system's equilibrium thermodynamic functions. 

The EA theory was followed by numerous theoretical 
studies that introduced new ideas, concepts and novel 
mathematical frameworks aimed at clarifying and un- 
derstanding the unconventional nature of the spin-glass 
phase transition, the low-temperature phase and the EA 
order parameter 0-0]. Importantly, similar to the EA 
theory, current approaches are based on the existence of 
a second distinct spin glass phase and a phase transi- 
tion of some sort, and discuss the associated order pa- 
rameters. Consequently, the predominant current view 
is that spin glass transition is a thermodynamic phase 
transition, even though unconventional [iHs, 0, S] ■ 
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Several outstanding problems remain with the phase 
transition approach. An essential problem rests with 
identifying the nature of the spin glass phase. This begs 
the important question of whether experimental data can 
be understood without making the assumption that a 
distinct spin glass phase exists. Another important is- 
sue is the conflict between theory and experiments [![. 
Experiments widely show strong dynamic effects. T g sig- 
nificantly increases with field frequency, or as observation 
time decreases [H, [jl-Eojl . This is inconsistent with a tran- 
sition between two equilibrium phases We note at 
this point that the same effect is seen in the liquid-glass 
(structural) transition [ill, [jij]. 

In addition to the susceptibility cusp, other proper- 
ties of spin glass systems have been studied, but have 
not lead to a conclusive evidence for a phase transition 
0, 0, 0] • For example, magnetic heat capacity shows no 
anomaly at T g , unlike in a phase transition 0] • On 
the other hand, the behavior of non-linear susceptibility, 
Xnh above T g is interpreted to be consistent with criti- 
cal scaling and is taken as evidence for the static phase 
transition at T g Q. This conclusion is shared by most 
authors @, although others are more cautious It 
is noted that scaling exists away from T g only, but fails 
close to T g where \ni flattens off, inconsistent with a 
phase transition picture [2, 0] • There are other problems 
with analyzing this and other similar scaling data: first, 
the choice of T g is arbitrary because T g is not fixed but 
depends on observation time, introducing arbitrariness in 
the scaling analysis. This applies not only to the interpre- 
tation of experiments, but also to extensive simulations 
of spin glasses 0. In particular, as larger times are ac- 
cessed in the simulations, T g decreases as expected from 
the experiments, altering critical exponents, lower criti- 
cal dimension and so on [8] , and thus suggesting that no 
firm values of these quantities can be asserted. Second, a 
phase transition is accompanied by the scaling behavior, 
but the opposite does not generally need to be the case. 
Consequently, the scaling behavior is seen in spin glass 
simulations where the phase transition is known to be 
absent [H, 0] , weakening the evidence for a phase transi- 
tion from the scaling data. Finally, the system is not in 
equilibrium at and below T g , so it is not known that the 
scaling data are related to a new phase 0, 0, 0] ■ We note 
that slow non-equilibrium relaxation effects could pre- 
cede a true phase transition as well (eg paramagnetic to 
ferromagnetic), however in this case the presence of the 
distinct low-temperature phase is known experimentally. 

The non-equilibrium state of the spin glass system gen- 
erally complicates the approach based on a phase tran- 
sition. In both experiments and simulations, relaxation 
time t continuously increases on lowering the tempera- 
ture [H, 0, 0] , as is the case in the structural glass tran- 
sition. At a certain temperature, r exceeds observation 
time, at which point the system falls out of equilibrium. 
This presents a problem for understanding experimental 



systems using statistical mechanics Q • This problem and 
other issues discussed above, including the increase of T g 
with field frequency, remain unresolved at present Q. 

Notably, profound non-equilibrium dynamic effects as 
well as thermodynamic anomalies are observed in both 
spin and structural glass transitions. Interestingly, the 
current view of the two phenomena is markedly differ- 
ent. The predominant opinion is that the structural glass 
transition is not a phase transition, and is a purely dy- 
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namic effect despite the heat capacity jump at T g 
The current view also holds that the spin glass transi- 
tion is a phase transition, the unconventionality of which 
comes from the dynamic effects [H, 0, 0, 0] . For exam- 
ple, Mydosh asks: "How can we get rid of the dynamical 
processes to determine the basic properties of the under- 
lying phase transition? At present an answer has not yet 
appeared" 

In this paper, we propose that instead of attempting to 
get rid of universal and profound dynamic effects in order 
to discuss the underlying phase transition, a more fruitful 
approach is to consider that the dynamic effects are at 
the heart of the spin glass transition. In this approach, 
we explore whether the main experimental results can 
be understood in an entirely dynamic picture. Contrary 
to the assumption currently held 0, we propose that 
paramagnetic and spin glass states are qualitatively the 
same. Consequently, we do not discuss a phase tran- 
sition of some sort, and therefore do not face complica- 
tions and controversies surrounding the nature of the spin 
glass phase and a phase transition [l|-0, 0] ■ We propose 
that the central experimental result, the susceptibility 
cusp, necessarily originates when a system stops relaxing 
at the experimental time scale. In this picture, the ob- 
served cusp marks the crossover that separates two time 
regimes, high-tcmperature relaxational regime t > r and 
low-tempcrature spin wave regime t < r, where t is the 
observation time and r is system relaxation time. Time- 
dependent effects, inconsistent with the phase transition 
approach, and the logarithmic increase of T g with field 
frequency in particular, originate as the immediate con- 
sequence of the proposed picture. Our proposal thus rep- 
resents a new outlook at the problem in which the exper- 
imental data can be interpreted in a way that is simpler 
and more physically transparent, an attractive possibil- 
ity in view of persisting difficulties and controversies in 
the field. 

In our discussion, we explore close parallels between 
the spin glass transition and the structural glass tran- 
sition. Importantly, we do so opposite to the historical 
trend when the insights from the existing spin glass the- 
ories were used to understand the structural glass tran- 
sition (see, e.g., Refs. d, 33"34|). These previous ap- 
proaches invoked the concepts from the spin glass transi- 
tion theories based on phase transitions and subsequently 
applied them to the structural glass transition. Conse- 
quently, the structural glass transition was related to a 
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phase transition of some sort. On the other hand, we 
start from a recent picture which explains the structural 
glass transition as a purely dynamic phenomenon, with- 
out any reference to a distinct solid glass phase or a phase 
transition 35]. We subsequently propose that the spin 
glass transition can equally be understood as a purely dy- 
namic phenomenon. The next section outlines the recent 
dynamic approach to structural glass transition, followed 
by our proposal that the same general idea can be used 
to understand the spin glass transition. 



STRUCTURAL GLASS TRANSITION 



Dynamic origin of the heat capacity jump at T g 



In the field of spin glass transition, the most commonly 



discussed effect is the cusp of susceptibilit; 
the area of structural glass transition 



at T a . In 



2fflXk, the re- 



sponse of a system is measured to the heat flow, and 
is quantified by heat capacity. If a liquid avoids crys- 
tallization on cooling, the heat capacity changes with a 
jump at the glass transition temperature T g . For var- 
ious liquids, ^ = 1.1 — 1.8 [25], where C l p and are 
constant-pressure liquid and glass heat capacities, respec- 
tively. The jump of C v provides one definition of T g at 
which the glass is said to form. Because above and be- 
low T g the system structure is equally disordered, the 
jump of C p immediately presents a problem that is at 
the heart of glass transition 0, (2||: how can the jump 
be understood if there is no distinct second phase? 

This problem remains unresolved, resulting in several 
important and interesting controversies 23|-29j]. One set 
of theories rationalize the jump of heat capacity by re- 
lating it to the existence of a phase transition and using 
the concepts from the phase transitions theory. This ap- 
proach has been convin cing ly criticized for a number of 
reasons (see, e.g., Ref. [Hf). Perhaps the most impor- 
tant reason is that it has not been possible to identify 
a distinct low-temperature phase, the glass phase. To 
get round this problem, several theories have put for- 
ward the proposals about the non-conventional mecha- 
nisms of the phase transition and non-trivial descriptions 
of the second phase, while retaining the idea of a phase 
transition in general [23| . Another set of glass transition 
theories maintain that glass transition phenomena at T g 
have purely dynamic origin, corresponding to the freez- 
ing of atomic jumps in a liquid at the experimental time 



scale [23j. The assumed absence of a phase transition 



and thermodynamic anomalies is supported by the wide 
experimental observation that the liquid and the glass 
structures at T g are nearly identical [23j. However, the 
challenge for these theories is to explain the origin of the 
heat capacity jump at T g as well as its large magnitude, 
which can be of the order of /cb per atom. 



We have recently proposed [3|| how to explain the 
jump of heat capacity in a purely dynamic picture, with- 
out asserting the existence of a distinct phase and a phase 
transition of some sort and, therefore, showed how to 
reconcile the above contradiction. We recall that glass 
transition temperature T g has two experimental defini- 
tions which give similar values of T g . In the calorimetry 
experiments, T g is the temperature at which the jump of 
C p is measured. In the experiments that measure liquid 
relaxation time r (e.g. dielectric relaxation experiments), 
T g is the temperature at which r exceeds the experimen- 
tal time t of typically 10 2 — 10 3 s [23|. We proposed that 
when r exceeds t, the jump of heat capacity at T g fol- 
lows as a necessary consequence due to the change of liq- 
uid elastic, vibrational and thermal properties including 
the bulk modulus and thermal expansion. Therefore, the 
jump of C p at T g is the result of the dynamic crossover 
rather than a phase transition. Consequently, there is no 
need to invoke a second glass phase, a phase transition 
of some sort and face associated problems. 

We first discuss why and how liquid thermal expansion 
coefficient a and bulk modulus B change at T g . Unlike in 
a solid, atoms in a liquid are not fixed, but rearrange in 
space. This gives the liquid its ability to flow. Each flow 
event is a jump of an atom from its surrounding atomic 
"cage" , accompanied by large-scale rearrangement of the 
cage atoms. We call this process a local relaxation event 
(LRE). A LRE lasts on the order of elementary (Debye) 
vibration period tq — 0.1 ps. Frenkel introduced liquid 
relaxation time r as the time between LREs in a liquid at 
one point in space 37[. Smaller and larger t correspond 



to smaller and larger liquid viscosities, respectively. Bas- 
ing on this property only, Frenkel concluded that at short 
times t < r, liquid response is the same as that of a solid, 
i.e. is purely elastic. On the other hand, for t > r, vis- 
cous liquid flow takes place, during which each LRE gives 
the additional, viscous, displacement. Hence, for t > r, 
liquid response to external perturbation (e.g. pressure) 
consists of both elastic and viscous response [37( . This 
picture provided the microscopic basis for the earlier phc- 
nomenological model of Maxwell [38|, who proposed to 
separate elastic and viscous response in his viscoelastic 
approach to liquids. 

Lets consider that a liquid is subject to pressure P. 
Pressure induces a certain number of LREs, which bring 
the liquid to equilibrium at new external conditions (P, 
T) after time r. As follows from the Maxwell-Frenkel 
approach, the change of liquid volume, v, is v — v e i + v r , 
where v e i and v r correspond to solid-like elastic deforma- 
tion and viscous relaxation process due to LREs, respec- 
tively. Lets define T g as the temperature at which r ex- 
ceeds the observation time t. This implies that LREs do 
not operate at T g during the time of observation. There- 
fore, v at T g is given by purely elastic displacement as in 
elastic solid. Then, P = Bi^^ and P = B g ^ above 
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and below T g , respectively, where Bi and B g are bulk 
moduli of the liquid and the glass, VJ° and V g are initial 
volumes of the liquid and the glass, and v g is the elastic 
deformation of the glass. Let AT be a narrow tempera- 
ture interval that separates the liquid from the glass so 
that t in the liquid, r u is t ; = r(T g + AT) and ^<1. 

Then, VJ° rs V g and v e i « v g . Combining the expressions 
for Bi and B g , we find: 



where ei = — is the ratio of relaxational and elastic 
response to pressure. 

The relationship between cti and a g can be calculated 
in a similar way, giving 35] : 

m = (e 2 + l)a g (2) 

where — — is the ratio of relaxational and elastic 
response to temperature. 

Eqs. (|H2p describe the relationships between B and a 
in the liquid and the glass due to the presence of LREs 
in the liquid above T g and their absence in the glass at 
and below T g , as long as T g is the temperature at which 
t < t. We note that B g and a g can be called unrelaxed, 
or non-equilibrium values, of bulk modulus and thermal 
expansion, respectively, because relaxational response v r 
in the glass decays during time r. 

We now relate the jump of constant-pressure heat ca- 
pacity, C p , to the corresponding changes of a and B 
at T g . We write C p = VTa 2 B + C v , where C v is the 
constant- volume heat capacity. We note that C v around 
T g , or at any temperature where r ^> To, is essentially due 
to the vibrational motion because the relative number of 
diffusing atoms is proportional to 21 , and is therefore 
negligibly small [35j. Using the Griineisen approxima- 
tion, we have shown that C v = 3N(1 + aT) [351 ] . 

As discussed above, a and B are different above T g 
when t > t and below T g when t < r (see Eqs. ([T|[2])). 
Therefore, the two different time regimes give different 
values of C p : 

t>r: C l = VTaf B x + 3N(l + a t T) (3) 



t<r: C^VTa 2 g B g + 3N{l + a g T) (4) 

where C l p and C| refer to liquid and glass, respectively. 

Eqs. (J3JU) relate the change of heat capacity to the 
changes of a and B due to the presence of relaxational 
response in the liquid and its absence in the glass. Using 
the experimental values of a g , ai, B g and Bi, we have 
shown that the calculated values of the liquid and glass 
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according to Eqs. (|3|4[) are in reasonable agreement 

with the experimental [351 . 

Importantly, the jump of heat capacity at T g in our 
theory takes place within the same single thermodynamic 
liquid phase, but below and above T g the liquid is charac- 
terized by different values of a and B due to the freezing 
of LREs at T g where the liquid falls out of equilibrium. 
Therefore, our theory is purely dynamic. In contrast to 
previous theories of glass transition E>ij , we do not dis- 
cuss transitions between distinct thermodynamic phases, 
even though it may be tempting to invoke phase transi- 
tions, conventional or unconventional, in order to explain 
the origin of heat capacity jump. Instead, the jump of 
heat capacity in this picture is due to the crossover be- 
tween two different regimes, relaxational regime above T g 
and elastic regime below T g . As we have shown earlier, 
the crossover may be sharp enough if small variations of 
temperature around T g give large changes of r, leading 
to the freezing of LREs in a narrow temperature range 
and giving the appearance of a phase transition [39} . 

As discussed in Ref. [35] in detail, our theory pre- 
dicts no difference of C p between the liquid and the glass 
at long times when both systems are in equilibrium. In 
this case, Bi = B g and ai = a g , giving Cp — Cp 1 . We 
note that reaching such an equilibrium state for com- 
mon glasses can take astronomical times and longer. 
For example, lets consider SiO*2 glass at room temper- 
ature T r =300 K. The activation energy barrier U can 
be assumed constant, because SiC"2 is a "strong" liq- 
uid [25] ]. Then, combining r(T g ) = TQCxp(U/T g ) and 

r(T r ) = roexp([//T r ), we write r(T r ) = r (^) 

Taking r =0.1 ps, T g 1500 K and r{T g ) = 10 3 s, we 
find r(T r ) = 10 67 s, approximately the fourth power of 
the age of the Universe. Although solid for any practi- 
cal purpose, Si02 at room temperature is an equilibrium 
liquid during times t > r(T r ). Consequently, it shows no 
jump of heat capacity on cooling from high temperature 
if t > r{T r ). 



Time-dependent effects 

The dynamic origin of the jump of C p at T g in our the- 
ory explains a well-known effect that T g logarithmically 
increases with the quench rate q = , where AT is the 
temperature interval above T g (see, e.g., Refs. pH I^L 
According to our theory, the jump of heat capacity at T g 
takes place when the observation time t crosses relaxation 
time r. This implies that r at which the jump of heat 
capacity takes place is r(T g ) = Combining this with 

r(T g ) = Toexp (jf-j (here U is nearly constant because 
r is approximately Arrhenius around T g [40| ]) gives 
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T„ 



U 



In 



AT 



— In q 



(5) 



Xc = -y~ (7) 
Curie behaviour of y is observed in spin glass systems 
at high temperature 

On temperature decrease, a spin-glass system can not 
find an ordered state due to frustration related to dis- 
order of various types [l[ . If C/ is the activation barrier 
for a spin rearrangement, r = roexp(^), where To is 
the elementary vibration period. U can be temperature- 
dependent, as discussed in the next section. Hence, r 
continuously increases on lowering the temperature until 
condition t > r is violated at a certain temperature, T g . 
When t — t at T g , spin LREs become frozen at the ex- 
perimental time scale. At this point, the system falls out 
of equilibrium and becomes non-ergodic. Consequently, 
Eqs. © and (J7J do not apply at and below T g . The re- 
maining excitations in the system are spin waves. Asso- 
ciated, in classical representation, with small-amplitude 



According to Eq. ([5]), T g increases with the logarithm 
of q, as widely observed in the experiments. Further- 
more, the dependence o f T a on q is consistent with the 
experimental results 2l|,|22j). This is the immediate and 
important result of our theory. 

We note that Eq. §5§ predicts no divergence of T g 
because the maximal physically possible quench rate is 
set by the minimal elementary time To (Debye vibration 
period) so that ^ is always larger than q in Eq. (JS|). 



SPIN GLASS TRANSITION 

Our main proposal is that, similar to the jump of B, 
a and C p in the structural glass transition, the suscep- 
tibility cusp in the spin glass transition is not related 
to the existence of a distinct spin glass phase and a 
phase transition of any sort. Instead, we propose that 
the cusp is the result of the crossover between the two 
regimes, the high-temperature relaxational regime and 
low-temperature spin wave elastic-like regime, when ob- 
servation time t exceeds system relaxation time t. In 
this picture, the crossover is a purely dynamic effect that 
takes place within the same single phase. 



Susceptibility above and below T g 

Similar to liquids, we define t as the time between 
two consecutive spin jumps (large-scale spin rearrange- 
ments in classical picture or transitions between states 
with different spin quantum numbers) that play the role 
of LREs in a liquid. When t > r at high temperature, 
spin LREs are operative, governing the dynamics of the 
system. In this regime, the system is in equilibrium that 
is equivalent to the equilibrium relaxational state in a liq- 
uid at T > T g . Consequently, the equilibrium statistical 
mechanics can be applied to spin LREs. In a two-level 
system of non-interacting spins, the free energy is 4l|: 



F 



-NT In 2 cosh 



fj,H 
~T~ 



(6) 



where /i is the spin magnetic moment, H is the applied 
field and fee = 1- The magnetic moment M = — ^jj — 
Nji tanh , which at high temperature becomes M — 



T ' 



spin displacements [41[ , spin waves carry thermal energy 
in an interacting system, ordered or disordered. In spin 
glasses, spin waves are analo gues of elastic waves prop- 
agating in glasses or liquids 37|, and can be similarly 



localized and damped at short wavelengths Q. At large 
wavelengths, spin waves in spin glasses have been pre- 
dicted theoretically, with the linear dispersion u = ck, 
where c is the speed of magnons 42[. This has sub- 



sequently been confirmed in theoretical studies 

HIT 



simulations |45| and experiments 



43|,|44|, 



spin waves are present above T„ as well 



We note that 
but this does 



not alter our main results because spin waves become 
increasingly damped at high temperature [42| . 

Unlike magnetic LREs that fall out of equilibrium at 
T g , the gas of magnons is in thermal equilibrium. Similar 
to the free energy of phonons 0] , the free energy of non- 
interacting magnons is: 



F = Ne + T ^2 In f 1 - exp 



flbJi 



(8) 



where are magnon frequencies and Neo is the zero- 
point energy. 

Magnetic moment and susceptibility due to spin waves, 
Xswi have been calculated for ordered magnetic struc- 
tures. For example, using the dependence of magnon 
frequencies in a ferromagnet on external field H, huik = 
2JSk 2 a 2 + qhbH, where J is the exchange parameter, S 
is spin and a is lattice constant in Eq. (JSJ, gives Bloch 
law: M oc -T§ 

On the other hand, Xsw has not been calculated for 



giving Curie law for susceptibility xc = 177 



spin glasses. From Eq. 
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where xo is the susceptibility due to zero-point vibra- 
tions. 

For spin glasses, field dependencies of u>i can be taken 



as loi 



ck and to t 



{ckf 



(if) ±3r, where w , 
and Wt are longitudinal and transverse waves, respec- 
tively 0. This gives {^f = (§-) 2 and *g 



4h 2 



for small fields (H — > 0) at which x is measured in spin 
glasses [l|. Next, the linear u) oc fe relationship for small 
fields implies the quadratic density of states, g(u>), as for 



phonons. g(co) 



6jy w 2 



where is Debye frequency 



and the normalization coefficient takes into account that 
2N transverse waves contribute to the sums in Eq. © 
because field derivatives of the longitudinal frequency are 
zero and hence do not contribute to \sw ■ Using the above 
field derivatives, substituting the sums in Eq. (0) with 
integrals with g{ui) and noting that integration can be 
extended to oo at low temperature due to the fast con- 
vergence of the integrals, we find: 



Xs 



Xo + 



TT 2 2 T 2 

4 " T£ 



(10) 



where Trj = huj is Debye temperature of spin waves. 

We find, therefore, that Xsw is quadratic with temper- 
ature and approaches constant value xo at T — > 0, consis- 
tent with experimental results (see, e.g., Refs (H [Til Il3l. 
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53j). Interestingly, the same behaviour (quadratic at 
low temperature and temperature-independent at T — > 
0) was found in the EA theory based on the description 
of the spin glass transition in terms of a phase transition 
and order parameters [gj]. 

We note that the approximations used above, includ- 
ing the limit of zero H and the dependence of uj on k 
and H, may affect the low-temperature behavior of x- 
However, it is clear that the behavior of Xsw below T g is 
qualitatively different from xc above T g due to a different 
mechanism of magnetic response. 



Dynamic crossover at T g 

As long as spin LREs are operative and in equilibrium 
when t > t, Eqs. © and (JT)) apply. This gives x = Xc 
in the long-time relaxational regime above T g : 



Regime (|11|) is equivalent to the high-temperature 
regime in the structural glass transition where LREs in a 
liquid are in equilibrium when t > r. On the other hand, 
when t < t in the short-time spin- wave regime, magnetic 
response is governed by Eqs. ((SJ and (flO)) . This gives 



X = Xs 



at and below T„: 



t <t : x = Xs 



(12) 



t > t : X = XC 



(11) 



Therefore, we find that temperature decrease results 
in the dynamic crossover between the two regimes, relax- 
ational (fTTj) and spin- wave regime (fT"2")) . This crossover is 
analogous to the crossover between relaxational and elas- 
tic regime in liquids, see Eqs. ([3]) and (jlj). We propose 
that the experimental susceptibility cusp is due to this 
crossover. We make two further observations supporting 
this proposal. 

First, we show that the crossover at T g can be fairly 
sharp, giving the appearance of a "cusp" characteristic 
of true phase transitions. We identified spin glass sys- 
tems for which the temperature dependence of r has been 
measured, and selected several representative systems in 
a wide range of T g between about 3 K and 300 K (see 
Table 1). In Refs. r has been fitted to the Vogel- 

Fulcher-Tammann (VFT) law: r = tq exp (^ T A To J . Us- 
ing the dependence t(T), we calculate what temperature 
increase above T g , AT g , is needed in order to reduce r 
by an order of magnitude: r(T g ) = 10t(T 5 + AT g ). AT g 
gives an estimation of the temperature range in which the 
predicted crossover between the high-temperature relax- 
ational regime above T g and low-temperature spin-wave 
regime at T g operates. The values of both Tf and AT/ 
are given in Table 1. We also list parameters A, To and 
T g used in the calculation. 

According to Table 1, AT g varies in a narrow tempera- 
ture interval relative to T g , from a fraction of K to several 
K in systems with small and large T g , respectively, giv- 
ing -=5* w 2% on average. Consequently, the crossover 
between the two regimes can be sharp and appear cusp- 
like. Further, lets define AT X as the approximate exper- 
imental temperature range around T g in which the high- 
temperature Curie-like behaviour of x crosses over to the 
low-temperature spin-glass dependence. AT X is shown in 
the last column of Table 1, and we observe that AT X is 
generally much larger than AT g . This means that in the 
temperature range where the experimental crossover of x 
operates, r increases by many orders of magnitude. This 
increase is more than sufficient to induce the proposed 
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TABLE I: T g , AT g and AT X for several spin-glass systems. 



System 




To 


T a 


ATg 


ATg 


AT X 




(K) (K) 


(K) 


(K) 




(K) 


Nao.3NiO2-l.3H2O [9] 


72 


1.9 


3.5 


0.1 


2.5 


2 


CdCri.9Ino.1S4 [10] 


114 


6.1 


10.6 


0.4 


4.2 


2 


U 2 CuSi 3 [11] 


57 


17 


19 


0.2 


0.9 


3 


CdCri.8lno.2S4 [10] 


132 


14.1 


21.1 


1.0 


4.7 


2 


CuMn (4.6 %) [12] 


59 


25.5 


27.5 


0.2 


0.6 


2 


LaMno.7Mgo.3O3 [13] 


100 


27 


30.5 


0.3 


1.0 


4 


FegsZry [14] 


300 


100 


105 


0.2 


1.9 


10 


Cuo.2Zno.8Fe 2 4 [15] 


85 


100 


110 


3.7 


3.4 


15 


Nio.7Zno.3Fe1.7Tio.3O4 


[16] 50 


218 


226 


4.7 


2.1 


75 


PMN [17J 


472 


291.5 


311.9 


2.3 


0.7 


10 



T(K) 



16 



20 



24 



FIG. 1: Experimental spin glass susceptibility [llj showing 
the crossover between high- and low-temperature behaviour. 



crossover between free-moving spins in the paramagnetic 
regime and frozen spins in the spin-wave regime. 

Therefore, Table 1 illustrates that (a) small temper- 
ature variations around T g are enough to give order- 
of-magnitude changes of r and induce the proposed 
crossover between regimes (fTTT) and (fL2"]) . and (b) the 
proposed mechanism can explain the experimental be- 
haviour of x because the predicted temperature range of 
the crossover is consistent with the range in which the 
experimental crossover (or, if appropriate, the cusp) of x 
develops. 

We observe at this point that the character of change 
of x around T g in many spin glass systems warrants the 
term "crossover" rather than the "cusp", as witnessed 
by smooth and gradual experimental behavior around T ? 
and large values of AT X in particular H 11 -IB, H 20, 
531 ] . Figure 1 illustrates this point. 

Second, if the observed susceptibility cusp is due to 
the crossover between two regimes, it is interesting to see 



that xc and Xsw attain similar values at the crossover 
temperatures, T c , that are typically observed in the ex- 
periment. Order-of- magnitude estimation of T c can be 
done by comparing Eq. ([7]) and Eq. (|10j) and ignor- 
ing system-dependent xo m Eq. (|10l) . a good approxi- 
mation at higher temperature when xo is small relative 
to the second temperature-dependent term. This gives 

T c = To (4t)* « 0.7T D . If typical T D of spin waves are 
in the range between tens and hundreds of K (e.g. 40-70 
K [H Esj]), T c falls in the range of experimental values 
of T g (see Table 1). 

In the proposed picture, the susceptibility cusp, or 
crossover, is related to the freezing of spin LREs at the 
experimental time scale, similar to the structural glass 
transition. Consequently, the experimental data can 
be understood as not a result of a phase transition of 
some sort between two distinct magnetic phases as dis- 
cussed previously (HH, 0] ■ Instead, the effects at T g take 
place within the same single thermodynamic paramag- 
netic phase, but the system response crosses over from re- 
laxational above T g to spin- wave below T g . In this sense, 
our theory is purely dynamic. In contrast to previous 
spin glass transition theories [IH! , we do not discuss 
thermodynamic phase transitions and order parameters, 
even though it may be tempting to invoke phase transi- 
tions, conventional or unconventional, in order to explain 
the cusp of x- Consequently, our approach involves fewer 
assumptions and does not require new concepts such as 
broken ergodicity, ultrametricity, replica symmetry, com- 
plexity of the energy landscape and so on in order to ex- 
plain spin glass transition [JQ. Further, in contrast to 
previous theories, our approach does not face the prob- 
lems related to the presence of profound relaxation effects 
in spin glass transition discussed below. 

We note that the time-dependent origin of the 
crossover of x implies that the crossover disappears when 
the spin glass system becomes an equilibrium paramag- 
net at t > r. This is analogous to the disappearance of 
the heat capacity jump in a liquid when t > r. However, 
reaching the equilibrium regime at low temperature re- 
quires times significantly exceeding the duration of the 
measurement (see Si02 example in Section 2). 



Time-dependent effects 



Spin glass systems universally sho w p rofound time- 
dependent relaxation effects [IHE 0,0, El El El 



16, 18 



Perhaps the most striking effect is the increase of T g 
with the logarithm of field frequency. Inconsistent with 
phase transition theories, these effects have remained a 
challenge for spin glass theories based on phase transi- 
tion and order parameters. To circumvent the problem, 
these theories added spin dynamics on top of the phase 
transition framework. This was done by, for example, 
incorporating relaxation effects into the spin-glass order 



8 



parameter (see, e.g., Ref. 5j|), introducing spin dynam- 
ics according to Glauber or Langevin equations of motion 
or postulating novel concepts such as hierarchy of relax- 
ation times, energy landscape or ultrametric structure 
of metastable states, with certain properties attached to 
each concept [3-0] • 

On the other hand, our dynamic approach does not 
require any extra steps to rationalize time-dependent ef- 
fects, and explains time-dependent effects in a simple and 
straightforward way. This has an added advantage of 
consistency, in that both the cusp (crossover) of x an d 
time-dependent effects are discussed on equal footing. 

As proposed, the susceptibility cusp is due to the dy- 
namic crossover between relaxational regime (jllj) . where 
t > t, and the spin- wave regime (|12p . where t < t. 
Therefore, the temperature at which the cusp takes place, 
T g , is defined from the condition t — t. Lets take r in the 

form of, for example, the VFT law: r = tq exp ( t ^ Tq J . 
Then, noting that observation time t is the inverse of field 
frequency i>, t = t gives 



T g = To + — 

in vq — In v 



(13) 



where vq — ^ is Debye frequency. 

We find, therefore, that the immediate consequence 
of the proposed theory is the increase of T g with In v. 
We further observe that Eq. (fT"3"|) has precisely the form 
observed in the experiments (see, e.g., Ref. [12j]). This 
is the important result of our theory. We note that the 
increase of T g with In v also follows if a simple Arrhenius 
activation form of r is assumed instead of the VFT law, 
giving T = in Eq. 

Eq. (flU)) is the analogue of Eq. (JSJ which describes 
the increase of T g with the quench rate in the structural 
glass transition. Similarly to the universality of Eq. (|5|) , 
the increase of T g with \nv according to Eq. (fT3"| is a 
universal phenomenon, widel y o bserved in spin glasses 
(see, e. g., Refs. [E [1, Ell [II 111, El, El^i^l ) ■ 



Non-linear susceptibility 



We comment on the behavior of the non-linear sus- 
ceptibility, Xni = 1 — 5^717 > w here \i is the linear sus- 
ceptibility in the limit H — )• [3|. Experimentally, 



logx„; oc -7loge, where e 



T-T a 



and 7 is close to 2 



[3, 0, ■ This is widely taken as the evidence for critical 
scaling and a phase transition [USB]- ^ tne same time, 
several important problems exist in interpreting this re- 
sult, including the arbitrariness in choosing T„ and the 
presence of non-equilibrium state around T g [2, • I m ~ 
portantly, the scaling of Xnl works only away from T g . On 
the other hand, log(xru) flattens off close to T g where the 



scaling fails, inconsistent with a phase transition picture 
[30- 

The above behavior of Xnl is qualitatively consistent 
with the proposed dynamic approach to spin glass tran- 
sition. Indeed, as discussed above (see Table 1), small 
temperature increases above T g result in the rapid de- 
crease of r. Therefore, the condition for the equilibrium 
paramagnetic regime, t > r, starts applying not far above 
T g . In this regime, M = Nfitanh^f, giving Xl = 
for small H . Using M and xi in Xnl and expanding tanh 
up to the third order in corresponding to the first 
non-linear term in Xnl , gives 



Xnl 



(14) 



Hence, log Xnl K — 21ogT. If written as a function 
of e = T 9 , where T g is an arbitrary temperature, 

log Xni = 1°S | ( i^") ~ 21og(e + 1). For temperature 
away from T g where e is sufficiently large, log Xnl is ap- 
proximately proportional to — 21oge, as in the experi- 
mental data. Closer to T g where e is small, logXnZ flat- 
tens off because log(e+ 1) tends to 0, consistent with the 
experiments [3,0. 

We note that Eq. (TTJ} and the above analysis are 
valid in the temperature range where t > r, whereas 
for smaller temperature and e they become approximate. 
For smaller e, Xnl becomes a non-equilibrium property 
that describes relaxation of induced magnetization in an 
experimental time window, and therefore depends on T 
or e as is seen experimentally. The detailed discussion of 
this effect is outside the scope of this paper. 



Relationship to spin Hamiltonians 

Existing spin glass transition theories typically start 
with a spin Hamiltonian of the type 0-0] : 



H 



4$> 



Si " Si 



(15) 



It is therefore instructive to restate the above idea with 
a reference to spin Hamiltonians. 

Using Hamiltonian (|15p . spin glass transition theories 
propose recipes to calculate thermodynamic functions 
and statistical averages in a spin glass system. Impor- 
tantly, the EA spin glass order parameter, q = (sp's- 2 " 1 ) 
or more complicated order parameters, are considered to 
be such statistical averages, and are calculated using the 
conventional statistical physics formalism [3-0, 0. Im- 
plicit in this procedure is the assumption that q, although 
originally defined as a dynamic correlator, also describes 
a distinct equilibrium spin glass phase. On the other 
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hand, we proposed that in order to explain experimental 
results, this assumption is unnecessary. Consequently, 
our approach enforces no order parameter and phase 
transition on the Hamiltonian. We suggested that the 
spin glass transition is a purely dynamic process where, 
as in Frenkel's theory of liquids, all correlations decay 
after time t > r 37jj . This approach implies q — at any 
finite temperature when t > r. 

Notably, our dynamic approach is based on the intro- 
duction of r from the outset. At the same time, r and the 
subsequent dynamic treatment could in principle be de- 
duced from the first-principles Hamiltonian (1151) . Indeed, 
the Hamiltonian specifies all interactions in the system 
that govern all of its dynamic properties such as spin 
LREs, their activation barriers and relaxation times as 
well as spin waves. However, understanding dynamic be- 
haviour of the system based on the description by the 
first-principles Hamiltonian such as (|15|) is a very com- 
plex task and arguably is not feasible for a real system, 
particularly for a disordered frustrated system such as 
the spin glass. In particular, this applies to the deter- 
mination of relaxation times which are governed by the 
activation barriers, and which in turn are governed by 
Jij in a complex way: with modern computers, it is pos- 
sible to calculate the energy landscape of an interacting 
system with up to a hundred particles only. 

In this sense, the dynamic approach discussed here is a 
reduction, in that r is introduced, but not derived from 
the general Hamiltonian. Similarly, Frenkel's approach 
to liquids 37| was a reduction from the first-principles 
description of a liquid by a general Hamiltonian such as 
(1151) , where atomic coordinates u are used instead of spin 
variables S: H — ^ ^ JijUiUj. This reduction surpassed 

ij 

the above complexity related to the first-principles treat- 
ment of the Hamiltonian, and was based on the empirical 
observation that liquids flow and that r, or viscosity, is 
the basic property of the flow. 

The introduction of r highlights another important 
difference between the dynamic description of spin-glass 
transition and previous thermodynamic approaches. The 
thermodynamic approach ignores the dynamics from the 
outset: it assumes that the system visits all of its states 
and therefore considers infinite waiting times. Implicit in 
this approach is the absence of r as the flow property of 
the system. 



SIMILARITIES BETWEEN THE STRUCTURAL 
AND SPIN GLASS TRANSITION: SUMMARY 

We have summarized the important features that are 
common to the structural and spin glass transition in 
Table 2. Points 1-9 in the Table refer to the similarities 
between liquids and magnetic systems, including effects 
at T g . Points 10-12 refer to the similarities of relaxation 



of both systems above T„, discussed in the next section. 



RELAXATION EFFECTS ABOVE T g 

As in structural glass transition, two main effects that 
set in above T g in spin glass systems are the VFT law 
mentioned above and stretched-exponential relaxation 
(SER). In liquids, slow relaxation sets in at high tem- 
perature T' c where r is on the order of picoseconds 55- 
157] . Above T' c , relaxation is exponential. It was assumed 
or postulated that SER is due to some sort of coopera- 
tivity of molecular relaxation in a liquid that sets in on 
temperature decrease, but the physical origin of cooper- 
ativity has remained unknown. 

We have recently proposed [40| that the physical ori- 
gin of the VFT law and SER in a liquid is the elastic 
interaction between LREs. Large atomic displacements 
due to the LRE distort the surrounding liquid and induce 
elastic waves. The wave frequency, lo, is on the order of 
Debye frequency because the wavelength is on the or- 
der of interatomic separations, implying that lo > 1/r in 
the range of r relevant for glass transition. As discussed 
by Frenkel [37], waves with frequency lo > 1/r propa- 
gate in a liquid as in a solid. The waves distort cages 
around other LRE centres in the liquid, and therefore af- 
fect their relaxation. Hence, we identified elastic waves 
as the physical mechanism of mediating interaction, and 
proposed that this interaction sets the cooperatiyity 
of relaxation whose origin was discussed [23M2O . 0, [29| 
but not understood from the physical point of view. The 
key question is the range of this interaction. As discussed 
in detail [4(|, this range is given by d c \: 



d c i = CT 



(16) 



where c is the speed of sound. 

The non-trivial point is that d c \ = ct increases with t. 
This is directly opposite to the commonly discussed decay 
of hydrodynamic waves, whose propagation range varies 
as 1/r. The difference is due to the solid-like regime of 
wave propagation, lot > 1, which is qualitatively different 
from the hydrodynamic regime, lot < 1. We called d c i 
liquid elasticity length because it defines the range over 
which two LREs interact with each other via induced 
elastic waves. Importantly, d c \ = ct increases on lowering 
the temperature because r increases. We proposed that 
this is the key to the origin of slow relaxation, the VFT 
law [U as well as SER pj . 

As such, the above mechanism is general enough to 
describe relaxation in any dynamic system where local 
entities relax and induce mediating waves (elastic, spin 
or other) that affect relaxation of other entities in the 
system. Hence, we propose that the same mechanism 
operates in spin glass systems. Here, a localized spin 
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TABLE II: Comparison of the main effects of structural and spin glass transition 



Property 


Structural liquids Spin glass systems 




and glasses 




1. Inability to find an ordered state at low temperature 


Yes 


Yes 


2. Relaxation is governed by the dynamics of activated local events 


Yes 


Yes 


3. Elastic (spin wave) response at short time t < t 


Yes 


Yes 


4. Relaxational response at long time t > r 


Yes 


Yes 


5. The crossover between equilibrium (ergodic) relaxational and non-equilibrium 






(non-ergodic) elastic regime at T g when t = r 


Yes 


Yes 


6. Anomalous behaviour of C p in liquids and \ m s P m systems at T g 


Yes 


Yes 


7. Logarithmic increase of T g with cooling rate and field frequency 


Yes 


Yes 


8. Distinct glass phase at and below T g 


No 


No 


9. Phase transition at T g 


No 


No 


Phenomena above T g : 






10. Mediating interactions and propagating waves 


Elastic 


Magnetic 


11. At high temperature, relaxation is exponential, non-cooperative and Arrhenius Yes 


Yes 


12. Crossover to slow non-exponential (stretched-exponential or logarithmic) 






and non- Arrhenius relaxation at low temperature 


Yes 


Yes 



rearrangement induces a spin wave that affects relaxation 
of other spins in the system. When d c i < d, where d is the 
characteristic distance between the neighbouring spins, 
spin LREs do not interact, resulting in exponential and 
Arrhenius relaxation as expected. However, when d c i > d 
on lowering the temperature, interaction between spin 
LREs sets in. Similar to liquids, we propose that this 
interaction is responsible for the VFT law and SER seen 
in spin glass systems 0, 0, @-Ht| ■ 

Eq. (fLo|) predicts that, similarly to liquids (lo] ]. two dy- 
namic crossovers in spin glass systems operate at n = — 
and T2 = — , where c is the speed of spin wave and L is 
system size. The first crossover is from exponential (non- 
cooperative) to non-exponential (cooperative) relaxation. 
The second crossover is from the VFT to a more Ar- 
rhenius relaxation when d c \ reaches system size L [icj . 
In contrast to structural glass transition, the dynamic 
crossovers in spin glass systems have not been experi- 
mentally studied in detail. Therefore, our prediction can 
be investigated in future high-temperature experiments 
in spin glass systems. We note that this picture predicts 
that r and T g increase with L when g? c i > L, provided 
that the thermalization length of the spin wave is larger 
than d c \ . Seen in supercooled liquids [59l. |60] | , this effect 
is also seen in spin glass systems [6ll465l | . Consistent with 
the proposed picture, it also supports the possibility of 
the second crossover at r 2 = -|. 

SUMMARY 

In this paper, we have proposed to take a new outlook 
at the spin glass transition problem, in view of persisting 



difficulties and controversies in the field. We proposed to 
interpret existing and future experimental data in a way 
that is simpler and more physically transparent. There 
is certainly more work to be done to clarify how the pro- 
posed ideas apply to the wealth of spin glass systems and 
phenomena. 

We have argued that the structural and spin glass tran- 
sition are similar, but in a different way than previously 
thought. Instead of asserting that the spin glass transi- 
tion is related to a phase transition and order parame- 
ters and rollingout this approach to the structural glass 
transition [30l434| . we proposed that the structural glass 
transition can be understood as an entirely dynamic phe- 
nomenon, and subsequently suggested that equally can 
the spin glass transition. 

We generally remark that a phase transition between 
two distinct equilibrium phases results in the change of 
system properties, a well-studied and understood topic 
in physics 0, [Bj]. However, there are many examples, in- 
cluding in structural and spin glass transition, where the 
system equally changes many of its important proper- 
ties at a certain temperature, yet the second phase and 
the order parameter are not apparent, the nature of the 
phase transition is not clear and the system shows pro- 
found dynamic effects. We have proposed that in this 
case, the change of system properties can be understood 
as an entirely dynamic phenomenon with no reference to 
a phase transition of some sort and order parameters. In 
this simple picture, system properties change at the dy- 
namic crossover when the system stops relaxing at the 
experimental time scale when t = r. 

Finally, we have explored closed similarities between 
relaxation properties of liquids and spin glass systems 
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above T g and proposed that slow relaxation effects are 
due to interactions between spin LREs via induced spin 
waves. This can stimulate future experiments to study 
high-temperature relaxation effects in spin glass systems. 

I am grateful to J. S. McCloy, V. V. Brazhkin, V. 
Heine, V. Ryzhov, and J. C. Phillips for discussions, and 
to EPSRC and SEPnet for support. 
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